function [dem_data, dem_param] =  readDem(dem_path, dem_par_path)

    % 读取DEM参数文件
    dem_param = gamma.readDemPar(dem_par_path, true);
    assert(isequal(dem_param.m_dem_proj, 'EQA'));
    assert(isequal(dem_param.m_ellipsoid, 'WGS 84'));
    assert(isequal(dem_param.m_datum, 'WGS 84'));

    % 读取DEM数据
    dem_data = gamma.readRealData(dem_path, dem_param.m_size_x, dem_param.m_size_y, dem_param.m_data_format);

end